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Abstract 

We distinguish three mechanisms of doping graphene. Density func- 
tional theory is used to show that electronegative molecule like F4- 
TCNQ and electropositive metals like K dope graphene p- and n-type 
respectively. These dopants are expected to lead to a decrease in carrier 
mobility arising from Coulomb scattering but without any hysteresis 
effects. Secondly, a novel doping mechanism is exhibited by Au which 
dopes bilayer graphene but not single layer. Thirdly, electrochemical 
doping is effected by redox reactions and can result in p-doping by 
humid atmospheres and n-doping by NH3 and toluene. 

1 Introduction 

The type and concentration of carriers in graphene can be controlled 
by adsorbates or substrates which transfer charge to graphene. p-type 
electronic doping results from a disparity between the work function 
of graphene and the electron affinity of the adsorbate. In the sim- 
plest case this is caused by a difference between the highest occupied 
molecular level (HOMO) and the lowest unoccupied molecular level 
(LUMO) of the adsorbate and graphene. This alignment as we shall 
see can depend on the substrate as we shall show that Au dopes bilayer 
graphene but is ineffective for single layer graphene. Electrochemical 
doping results from redox reactions near the graphene surface. The 
timescales for the two types of doping is different since electronic dop- 
ing occurs spontaneously while electrochemical doping requires longer 
times to overcome reaction and diffusion barriers. Hysteresis effects 
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for graphene based field effect devices can also be expected in electro- 
chemical doping since the concentration of the reactants can change at 
a rate slower than the change in gate voltages. 

Electropositive elements like metals are good candidates for elec- 
tronic n-type doping of graphene. Molecular beam deposition of tran- 
sition metal clusters, Ti, Fe and Pt, has shown that these metals effect 
n-doping of graphene whereas bulk deposition of Pt dopes graphene p- 
type[T]. Electronic doping on graphene was also demonstrated with K 
atoms deposited at low temperatures (20K) on graphene in ultrahigh 
vacuum (UHV)[2 . Synchrotron-based high resolution photoemission 
spectroscopy (PES) showed that graphene can be made p-type by the 
deposition of tetrafluoro-tetracyanoquinodimethane (F4-TCNQ). Be- 
cause of the molecular separation of the adsobate and graphene, all 
these dopants are expected to lead to a decrease in carrier mobility 
arising from Coulomb scattering without any hysteresis effects. 

Exposure of graphene to a humid atmosphere or gaseous nitrogen 
dioxide leads to p-type behaviour, while ammonia leads to n-type be- 
haviour [3] . A great surprise is that it was found, in particular experi- 
mental conditions, gaseous toluene leads to n-doping of graphene|l], in 
contradiction with theoretical calculations which predict no transfer of 
charge between graphene and the toluene molecule [12]. In this paper 
density functional calculations are used to investigate the transfer of 
charge between graphene and adsorbates. A electrochemical model for 
doping of graphene is also presented. 

2 Method 

The electronic ground state of graphene with different adsorbates was 
calculated using spin-optimised density functional theory as imple- 
mented in the AIMPRO codefSj [6]. The exchange correlation po- 
tential was described using the local density approximation (LDA). 
The core electron were treated using the Hartwigsen Goedecker Hut- 
ter (HGH) pseudopotencials[7]. The orbitals of the valence electrons 
consist of independent s—,p—,d— like Gaussian functions centred on 
atoms [5]. The electronic levels were filled using Fermi-Dirac statis- 
tics with kBT=0.01 eV and a metallic filling. The Brillouin zone 
was sampled with a grid of 8x8x1 k-points within the Monkhorst- 
Pack scheme[9] . Charge densities were Fourier-transformed using plane 
waves with an energy cut-off of 200 Ha. The spin-populations of the 
cell were optimised by starting from random spin distributions. We 
used two different graphene unit cells, 4x4x1 and 6x6x1, enclosing 
32 and 72 carbon atoms, respectively. During the relaxation all the 
atoms were allowed to move to their equilibrium positions. The in- 
fluence of the cell size was also studied with larger supercells and the 
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band structure obtained and the consequent charge transfer were very 
similar. 

3 Electronic doping 

The geometry of the relaxed F4-TCNQ molecule was found to be in 
good agreement with experimental results for its precursor TCNQ [10] . 
The electron affinity, the difference between the LUMO and the vac- 
uum level, was found to be 5.25 eV in agreement with an experimental 
measurement of 5.2 eV[TT]. This deep LUMO level when compared 
with the work function of graphene ~ 4.5 suggest that F4-TCNQ is a 
good candidate to act as an electronic acceptor on graphene. 

The F4-TCNQ molecule was placed on top of a 72 carbon atoms 
graphene supercell, 6x6x1, and during the relaxation all the atoms 
were allowed to move. In the minimum energy configuration the molecule 
remain parallel to the graphene layer with a interplanar separation of 
3.1 A. The binding energy was found to be 1.26 eV which is probably 
overestimated as a result of the LDA approximation. 

Figure[l]a) and b) shows respectively the electronic band structures 
of pristine graphene and graphene with an adsorbed molecule. Because 
of band folding, the Dirac point, where the tt and it* bonding and anti- 
bonding levels of graphene meet, occurs at F instead of the K point 
in the 6x6x1 graphene supercell used here. Comparison of the two 
electronic band structures shows that when the molecule is placed on 
top of graphene, the Fermi level is shifted to lie helow the Dirac point, 
and a new filled flat band appears which is related to the adsorbate. 
This suggests transfer of negative charge from graphene to F4-TCNQ 
This is supported by an analysis of the wave function shown in Figure [2] 
a) of the level marked A in Figure [l]b). This shows that the new filled 
band is strongly localised on the F4-TCNQ molecule. In addition, the 
empty level marked B, which was occupied in graphene, is delocalised 
over the graphene layer. This confirms that charge transfer between 
graphene and F4-TCNQ occurs. The position of the Fermi- level indi- 
cates about 0.3 electrons are transferred from graphene to a molecule 
of F4-TCNQ. 

To study the doping properties of metals, a single atom of Au and K 
respectiverly were placed on graphene in the three obvious adsorption 
sites: above the centre of the graphene hexagonal ring, above a C-C 
bond and above the carbon atom. In the case of Au the minimum 
energy configuration was found to be 2.27 A on top of a carbon atom 
whereas K sits 2.42 A directly above the middle of a hexagon optimising 
the number of K-C bonds. Thus the K atom is strongly bound to 
graphene with a binding energy estimated to be 1.51 eV in contrast 
with the 0.65 eV estimated for a Au atom. This is in part due to their 
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Figure 1: (Colour online) Band structure (eV) of a) pristine graphene and 
b) F4-TCNQ on top of graphene plotted in the vicinity of the Fermi energy 
along the high symmetry branches of the graphene Brillouin-zone. The 
Fermi level is set to zero. Full lines (red) denote occupied states while 
dashed lines (blie) show empty level4. The Fermi level is placed at zero. 
The curves shows the unoccupied F4-TCNQ levels around eV become 
occupied near the Dirac point indicating charge transfer. Colour online 



a) 




b) 

Figure 2: (Colour online) Plot of the wavefunction of the a) HOMO level 
at r marked A in figure [2]3) shows strong localisation on F4-TCNQ and the 
b) LUMO at r marked as B in figure [2^) shows strong delocalisation of a 
pi-bonding orbital over graphene but avoids F4-TCNQ. Red and blue lobes 
are of equal amplitude and opposite sign. 
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Figure 3: (Colour online) Spin-averaged electronic band structure (eV) of 
K on top of graphene in the vicinity of the Fermi energy. The Fermi level 
is set to zero. Full lines denote occupied states while dashed lines show 
empty levels. The bands around B, unoccupied for pristine graphene, are 
now occupied. 



different doping properties. Spin-optimization led to a spin-less ground 
state for K and a spin-1/2 state for Au. 

Figure [3] shows the electronic band structure of K on graphene. 
The Fermi level lies above the Dirac point which indicates doping of 
graphene with electrons. The wavefunction, shown in Figure HI a), of 
the empty flat level marked A in Figure |3] is strongly localised on the K 
atom and is derived from the 4s-level of K. The wavefunction (Figure 
|4]b)) of the filled level above the Dirac point, marked B is delocalised 
over the graphene layer. This confirms that K dopes graphene n-type 
and we estimate about one electron per K atom is transferred. 

We now turn to doping with Au. Figure [5] a) and b) show the ma- 
jority and minority spin electronic band structures for Au on graphene. 
Although there are changes to the dispersion of the bands in the vicin- 
ity of the Dirac point there is no significant change in the position of 
the Fermi level with respect to the Dirac point compared with pris- 
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Figure 4: (Colour online) Plot of the wavefunctions of the majority spin 
electronic levels of a K atom on top of graphene at points marked a) A and 
b) B in figure [31 The wavefunction of the level marked A is localised on the 
K atom while the wavefunction of the highest occupied level, marked B, is 
delocalised over the graphene layer. 
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tine graphene. This implies that there is no significant charge transfer 
between Au and graphene. 

To study the effect of additional graphene layers, a single Au atom 
was placed between the AB stacked graphene bilayer. The relaxed 
ground state was spin-less and the Au atom occupied a site nearly 
midway between the two graphene layers and 2.14 A directly above 
a carbon atom in the lower layer and directly below the middle of 
the hexagon of the top layer. The separation between the graphene 
layers increased to approximately 4.5 A. The graphene sheets do not 
appear to appreciably distort. Figure [5] a) and b) show the electronic 
band strucuture of a pristine graphene bilayer and a graphene bilayer 
intercalated with Au atoms. As for single layer graphene, the highest 
occupied and lowest empty bands touch at the Dirac point but the 
linear dispersion is lost. However, when the Au atom is placed between 
the graphene layers, the Fermi level is shifted above the Dirac point. 
An analyis of the wavefunction shows the flat level far from the K point, 
marked as A, results from the 6s Au level. When the 6s level crosses 
the TT* bands of graphene, they hybridise and charge is transferred to 
graphene effecting n-type doping. The level marked as B in Figure [SJ 
is delocalised over graphene and this level is now occupied, and unlike 
the case when the second graphene sheet was absent. We suggest 
that the effect of a second graphene layer is to compress the Au 6s- 
wavefunction with a consequent upward shift of the level. The same 
effect was observed when ther second graphene layer was replaced by 
a toluene molecule [12] . 

4 Electrochemical doping 

Another type of doping involves electrochemical redox reactions of 
graphene with water and adsorbates. It has been reported that ex- 
posure of graphene to a humid atmosphere or NO2 results in p-doping 
whereas exposure to gaseous NH3 [3] and toluene results in n-doping 
[4]. It is important to realise that other carbon based materials can 
exhibit similar effects. Diamond, whose surface has been treated with 
hydrogen which reduces considerable its work function, can also be 
doped p-type by exposure to a humid atmosphere [13]. The effect is 
suppressed by NII3 and enhanced by NO2 [15 • Similar effects are 
seen in carbon nanotubes [HI [16] . These types of doping cannot be 
explained by an electronic mechanism involving direct transfer of elec- 
trons from say toluene to graphene as the ionisation energy of toluene 
is ~ 6 eV [17] and greatly exceeds the work function of graphene. The 
transfer doping model has similar difhculties in the case of NH3 . How- 
ever, atmospheric doping has been suggested as being due to the direct 
transfer of charge between graphene and defects in layers of ice on top 
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Figure 5: (Colour online) Spin-polarised band structures (eV) of Au in the 
vicinity of the Fermi energy, (a) Majority spin band structure, and (b) 
minority spin band structure. The Fermi level is set to zero. Full lines 
denote occupied states while dashed lines show empty levels. 
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Figure 6: (Colour online) Electronic band structure (eV) for a Au atom 
intercalated in bilayer graphene in the vicinity of the Fermi energy. Full 
lines denote occupied states while dashed lines show empty levels. The 
Fermi level is placed at zero. Note tl^i^^band marked B which is unoccupied 
in bilayer graphene is now partially occupied showing doping has occurred. 



of the Si02 substrate [18]. Nevertheless, this mechanism could not ex- 
plain the surface doping of diamond as Si02 is not present. Instead we 
argue here that the doping effects due to the atmosphere, NO2 or N2O4, 
NH3 and toluene on graphene that lead to little change, or an increase, 
in carrier mobility are all examples of electrochemical redox reactions 
occurring in aqueous layers in contact with graphene. A similar theory 
will also be applicable to nanotubes and already has been invoked for 
diamond [19| . The similar properties of diamond films, nanotubes and 
graphene arise from their similar work functions although the large 
band gap of diamond excludes the possibility of n-type doping. 

The assumption that a water layer is in contact with graphene is 
a non-trivial one as theory has indicated that diamond and graphene 
are hydrophobic. We are inclined to the view that for graphene and 
nanotubes, likely places for water deposits are at the interfaces between 
graphene and Si02. Kelvin probe microscopy and X-ray spectroscopy 
have shown four or five water layers exist on the surface of Si02 films 
grown on Si [20] and there have been FTIR observations of overtones of 
vibrations of water molecules [21] It is also widely believed that water 
can exists in voids in a-Si02 [H] and would be resistant to thermal 
anneals if the voids are small. For diamond, we suppose water is found 
at places where the H termination has been replaced by electrophilic 
OH or oxygen groups. However, further work is required on this point. 

Redox reactions are reactions occurring in aqueous solution involv- 
ing changes in the charge states of the participants. The net charge 
change can be transferred to a electrode which in our case is graphene. 
Whether such reactions occur require that the change in the total Gibbs 
free energy is negative (described as a spontaneous reaction) and that 
the barriers are sufficiently small that the reactions can occur at room 
temperature. Tables of AG for a reaction in which unit charge is re- 
moved to or added from vacuum are available [23j but barrier heights 
are unknown and it must be assumed that they are sufficiently low. 
The total Gibbs free energy change is then AG -I- W for p-doping or 
AG — W for n-doping. Here W is the work function of graphene. 

It has been argued that the p-doping of diamond is due to the 
electrochemical reduction of O2 in the presence of water as in 

02+2H20 + 4:e^40H-. 

The reaction is spontaneous with the electron being transferred from 
the valence band of diamond and the formation of 0H~ groups jl9) . 
A current carrying hole is created near the diamond surface. We now 
show that an identical reaction is spontaneous for graphene. We con- 
sider the process occurs by removal of a electron from graphene to 
infinity followed by trapping of the electron by O2 and the subsequent 
reaction with water. The extraction of the electron from graphene re- 
quires the expenditure of energy equal to its work function W. This 
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has been measured for graphite in air to be 4.5 eV [24,. This value 
is close to 4.9 eV found for hydrogenated diamond [IS] and we expect 
the same work function for graphene j26| . The similarlities of the work 
function of diamond, graphite and nanotubes explains why the same 
redox theory pertains to all three materials. Now AG is the free en- 
ergy for the molecular reaction with the charge brought brought from 
vacuum and can be written as — ei? — 4.4 eV where E is the electrode 
potential relative to a hydrogen electrode. The reaction is assumed 
to involve molar concentrations of reactants and products. Standard 
electrode tables give E to be 0.42 V [53]. Thus the total free energy 
change for a electron to be taken from graphene and to reduce the 
02/water couple is -0.42-4.4 +W or -0.3 eV. Thus the reaction where 
electrons are extracted from graphene by oxygen and water is spon- 
taneous. The above assumes that that molar concentration of 0H~ 
and O2 are present but if this is not the case, then AG depends on 
the concentrations of 0H~ and O2 with the reaction being inhibited 
in strongly basic conditions. Nernst's equation gives AG to be 

AG = -4.8 - 0.059(14 - pH + 0.251ogio(po2))- 

Here P02 is the partial pressure of oxygen in atm. Taking this to be 0.21 
atm gives AG to vary from -4.8 eV for pH = 14 to -5.7 eV for pH = 
1. Acidic solutions will promote the reaction as the concentration of 
0H~ is so low. We note that 0H~ gradually builds up as doping 
increases and might be expected to lead to reduced carrier mobility 
as explained above. That this does not happen can be explained if 
0H~ is either very effectively screened by surrounding Si02 groups 
or is mobile through 8102- Previous calculations have shown that the 
latter is likely: the binding of OH^ to the network is only 0.3 eV 
and thus diffusion through Si02 could occur at room temperature [32] . 
Moreover, the theory explains the hysteresis effect [57] if the chemical 
reaction rates, or diffusion rate for water, oxygen and 0H~, are slower 
than the rate of change of gate voltage. For positive gate voltages, 
the reaction product 0H~ drifts away from graphene and the lower 
hydroxyl concentration there encourages further dissociation of water 
and hence greater increases in hole densities in graphene as indeed 
observed. [27] . 

The atmospheric doping effect disappears following vacuum anneal- 
ing around 200°C for both graphene and diamond [13]. We suppose 
this is because of the loss of oxygen and not to the loss of water. This is 
because we require a water layer to explain doping of annealed samples 
by gaseous NO2, NH3 and toluene. 

It might be thought that the effects of NO2 and NH3 can now be 
explained using the same redox couple but taking into account the 
changed pH. However this need not be the case and it is possible that 
different redox reactions occur especially in anaerobic conditions. A 
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possible reaction involves the reduction of N2O4 as in iV204 + 2e = 
2iVO- with a electrode potential of £; = .87 V [13]. Thus AG would 
be -5.3 eV and certainly could account for p-doping of both graphene 
and diamond. 

For the oxidation of NH3 or NH4OH we have, 2NH4OH ^ N2 + 
2H2O + 6H+ + 6e with E = -0.09 V and yielding AG = 4.5 eV. 
The electron produced by this reaction is transferred to graphene and 
the change in the total free energy is AG — W. Thus the oxidation 
reaction leads to the supply of electrons to graphene provided that the 
work function is greater than 4.5 eV which is just possible. However, 
we have ignored the basic conditions present as the reaction is inhibited 
in acidic conditions. Nernst's equation gives AG ranging from 4.5 eV 
{pH = 1) to 3.7 eV {pH — 14). Thus the energy cost in removing 
an electron from aqueous NH3 is 3.7 eV and is easily recovered by its 
transfer to graphene. The reaction occurs because of the scarcity of 
H"'" and a very high N2 bond energy. It is noteworthy that this reaction 
is used industrially to remove ammonia during waste water treatment. 

Finally we turn to toluene which has been found to result in n-type 
doping, a pronounced hysteresis effect and enhanced carrier mobil- 
ity [1]. One possibility is that this involves the redox couple where 
toluene is oxidised to benzyl alcohol. Such a reaction has been re- 
ported previously [55]. From the tables of free energies [33], AG for 
Toluene +20H- = Benzyl alcohol -I- H2O -I- 2e is 3.95 eV. Thus the 
electrochemical oxidation of toluene to benzyl alcohol is spontaneous 
as the work functions exceeds 3.95 eV and the liberated electron will be 
trapped by graphene achieving n-type doping. For negative bias con- 
dition, 0H~ will drift through the Si02 substrate towards graphene, 
removing a source of scattering centres and leading to an increase in 
mobility 

The reaction depends however on the concentration of OH^ in the 
Si02 which is very low. Another possibility is the anaerobic oxida- 
tion of toluene by water. This reaction is known to occur naturally 
at room temperature where it is catalysed by microorganisms. The 
standard electrode potential for Toluene -I- 21 H2O 7 IICO3" +43 
H"*" -t- 36e is 0.27 V [29]. In nature, the charge exchange takes place 
with the Fe^"'"/Fe^"'" couple [29]. Taking into account the standard elec- 
trode potential of Fc^+/Fe'^"'" is 0.77 V, this makes the total free energy 
change for Toluene + 21 H2O + 36 Fe^+ 7 HCO^ -h43 H+ +36 Fe'^+ 
to be -1.04 cV per electron. However, the reaction remains exothermic 
if the Fe^+/Fe'^+ couple is replaced by graphene, as then AG — is - 
0.37 eV. It is worth emphasising that redox reactions are slow and thus 
graphene needs to be exposed to toluene for long periods ^ 1 hour for 
a doping affect to be seen. Thus the anaerobic conversion of toluene 
into CO2 in the presence of graphene is favoured. To account for a 
mobility increase, we suppose that mobile H''" ions neutralise charged 
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scattering centres present in the undoped graphene. 

5 Conclusions 

In summary, we conclude there are two types of dopants for graphene. 
The first which can be called electronic doping occurs when there is a 
direct exchange of electrons with an adsorbate and graphene. A good 
example is K. Such doping occurs promptly and leads to a reduction in 
carrier mobility and there should be no hysteresis effects. The second, 
involving electrochemical doping occurs by redox reactions involving 
water at the interface. This can lead to an increase in carrier mobility 
but requires appreciable times to occur. This leads to hysterisis effects. 
The assumption that one or both of the charged products 0H~ or H+ 
is mobile and responds to the field due to the gate voltage could explain 
increases in carrier mobility. 
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